Mode-mixing quantum gates and entanglement without particle creation in 

periodically accelerated cavities 
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We show that mode-mixing quantum gates can be produced by non-uniform relativistic accel- 
eration. Periodic motion in cavities exhibits a series of resonant conditions producing entangling 
quantum gates between different frequency modes. The resonant condition associated with particle 
creation is the main feature of the dynamical Casimir effect which has been recently demonstrated 
superconducting circuits. We show that a second resonance, which has attracted less attention 
since it implies negligible particle production, produces a beam splitting quantum gate leading to 
a resonant enhancement of entanglement which can be used as the first evidence of the dynamical 
Casimir effect in mechanical oscillators. We propose a desktop experiment where the frequencies 
associated with this second resonance can be produced mechanically. 
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Introduction. — In relativistic quantum field theory, 
the particle content of a quantum state is affected by 
the evolution of the spacctime, including the motion of 
any boundaries. Further, the very notion of a "particle" 
depends on the motion of an observer. In flat space- 
time, celebrated examples are the thermality seen in 
Minkowski vacuum by uniformly accelerated observers, 
known as the Unruh effect P, 0] , and the creation of par- 
ticles by moving boundaries, known as the dynamical (or 
non-stationary) Casimir effect (DCE) 0, J. In curved 
spacetime, a celebrated example is the Hawking radia- 
tion emitted by black holes The DCE is related to a 
fundamental prediction by Fulling and Davies that a non- 
uniformly accelerated mirror will excite photons out of 
the vacuum [|| . It was later realised that this effect may 
be significantly enhanced if, instead of a simple mirror, 
a cavity is used in which one or both of the mirrors are 
in motion 0] • The simplest situation in which to observe 
the DCE is that of a cavity oscillating sinusoidally with 
frequency u r . The DCE is predicted to exhibit a funda- 
mental resonance condition for the production of quan- 
tum entangled photon-pairs, 0J r = u>\ + wj, where idxp 
are the two entangled photon frequencies [7| . The actual 
number of photons predicted for a mechanically oscillat- 
ing cavity is strongly limited (~ 10 -9 photons/second) by 
the maximum achievable oj r . For this reason a number of 
alternative systems that also exhibit a periodically vary- 
ing boundary of some kind have been proposed with the 
aim of enhancing the DCE. Examples are superconduct- 
ing SQUID mirrors, Bose-Einstein condensates (produc- 
ing phonon pairs) and cavities controlled using nonlinear 
optics [!, 0, IS ■ Notwithstanding recent breakthroughs, 
the DCE remains an extremely difficult effect to observe 



and study experimentally. 

In this letter we consider the general case of a rigid 
cavity undergoing an arbitrary (mechanically induced) 
acceleration. In the specific cases of a linear sinusoidal 
or a uniform circular motion, we show that the mode 
mixing resonance condition, u! r = \wi — W2I [III, can be 
brought significantly below the DCE photon generation 
resonance condition, to apparently experimentally acces- 
sible frequencies for which no new photons are generated. 
We show how this low-frequency resonance leads to the 
generation of entanglement between existing and previ- 
ously non-entangled cavity modes. The oscillating cavity 
can be shown to behave like a generalised beam-splitter, 
thus performing an essential quantum gate functionality. 
We then discuss the possibility of performing actual ex- 
periments with mechanically oscillating optical cavities. 

(1+1). — We first consider the simplified case of a 
cavity in (1 + l)-dimcnsional Minkowski spacetime. The 
cavity is assumed mechanically rigid, maintaintaining 
constant length L in its instantaneous rest frame. The 
proper acceleration at the centre of the cavity is denoted 
by <z(r), where r is the proper time. To maintain rigidity, 
the acceleration must be bounded by \a(T)\L/c 2 < 2 [111 ]. 
From now on we set c = H = 1. 

The cavity contains a real scalar field <j) of mass /io > 0, 
with Dirichlct boundary conditions. When the cavity 
is inertial, the field has a standard orthonormal basis 
of positive frequency mode functions, with the nondc- 
negenerate angular frequencies ui n = %/ M 2 + 7r 2 n 2 / L, 
11 = 1,2,..., where M = ^i L. We denote by a and 
/3 the Bogoliubov coefficient matrices between the stan- 
dard bases in an inertial initial region and a final in- 
ertial region, in the conventions of 11-13]. We write 
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U = ( ^ , so that the composition of Bogoliubov 
transformations is matrix multiplication of the corre- 
sponding [/-matrices, and the Bogoliubov identities (l3j 

read(i_° 1 )=f/(J-°i)C/ t - 

When the acceleration between the initial and fi- 
nal inertial regions is uniform and lasts for proper 
time f, we have Uh(f) = K^ 1 Zh(f)Kh, where 

o"h o/3 h \ ?t=.\ _ ( Z h{r) 



o z*(t))> ^( t ) 



diag(e 



Q, n (h) arc the angular fre- 



quencies during the acceleration, Q a^ and 0/ 8 h are the 
Bogoliubov coefficient matrices from the initial inertial 
segment to the uniformly accelerated segment, and the 
acceleration has been encoded in the dimensionless pa- 
rameter h = aL. 

For accelerations that may vary arbitrarily between 
the initial time To and final time r, U(t,tq) is given 
by the limit of U hN (f N )U tlN _ 1 (fjv-i) ■ ■ ■ U h2 (f 2 )U hl (n) 

as N — > oo, such that r — to = X^fcLi ^fc i s nx cd and 
each ffc — > 0. As an infinitesimal increase in r amounts 
to multiplying U(t,tq) from the left by Uh{f) with in- 
finitesimal f, U(t,tq) satisfies the differential equation 



U(t,t ) = iK h }Sl h ( T )K h ( T )U(T) 



(1) 



where ^Mt 



(t) 



'h(T) 




diag (Oi (/i(r)), f22(/i(T)), ••• ), and the overdot de- 
notes derivative with respect to r. The solution 
is 



u ( T fi T o) = Texp ( i / K h{T) VL h ( T )K h{T ) dr 



(2) 



where T denotes the time-ordered exponential. 

To summarise: the Bogoliubov transformation be- 
tween the inertial initial segment ending at proper time To 
and the final inertial segment starting at proper time tj is 
given by ([2]). h(r) may vary arbitrarily for tq < r < ry: 
in particular, no small acceleration approximation has 
been made. For piecewise constant h(r), @ reduces to 
a product of the matrices Uh{fk) from each constant h 
segment. 

A direct consequence of ^ is that the Bogoliubov co- 
efficients evolve by pure phases over any time interval in 
which h is constant. Particles in the cavity are hence cre- 
ated by changes in the acceleration, not by acceleration 
itself, as can be argued on general adiabaticity grounds 
[Til Hif . The cavity is in this respect similar to a sin- 
gle accelerating mirror, which excites photons from the 
vacuum only when its acceleration is non- uniform Q. 

(1+1): small accelerations. — At small accelera- 
tions, Q n (h), a a h and ft h have the expansions 

n n /u n = 1 + 0{h 2 ) , 71 = 1,2,..., (3a) 
a h = l + ha + 0(h 2 ), ft h = hp + 0(h 2 ), (3b) 



where the elementary, M-dependent expressions for a mn 
and ft mn are given in |ll| . a is antisymmetric and ft is 
symmetric, and both vanish when m — n is even. 
We seek U(t,tq) in the form 



a = e l " {T - To) 
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iw(r— To) 



-A + o(h 2 )), 
o(h 2 ), 



(4a) 
(4b) 



where u) = diag(u;x, w^, ■ ■ ■ ) and A and B arc of first 
order in h. Using ([!]) and ([3]), we find 



A mn = ih(T) (LU rn - Un ) e -^-^)(r-r ) 



B mn = ih{r) (oj m + cj n ) i 



-i(a; m +a; n )(T— To) 



, (5a) 

ftran ■ (5b) 



Once h(r) is specified, A mn and B mn may be found from 
([5]) by quadratures. 

Three observations are immediate from ([5]). 

First, A mn and B mn with even m—n remain vanishing, 
regardless of the acceleration. 

Second, in the special case of piecewise constant h, 
the changes in the magnitudes of A mn and B mn come 
entirely from the discontinuous jumps in h, as observed 
in (llj . However, if the jumps are replaced by a con- 
tinuous evolution that is slow compared with the os- 
cillating factors in ([5]), the changes in the magnitudes 
of A mn and B mn are significantly smaller. The as- 
sumption of rapidly-changing acceleration is hence essen- 
tial for the non-periodic travel scenario results obtained 
in [ll , 16 , 17 1 . We emphasise that no such rapid changes 
are involved in the experimental scenario that will be 
considered below. 

Third, if h is sinusoidal and its angular frequency 
oj r equals the angular frequency of an oscillating factor 
in ([5]) , the corresponding Bogoliubov coefficient will grow 
linearly in time. These resonance conditions read 



for A mri with odd m — 



n: ui r = \U>m — w 



m | 7 



(6a) 



for B mn with odd m — n: u r = uj m + u n . (6b) 



The particle creation resonance (I6bj 
DCE literature @, |, 0, M, 



is well known in the 
. The mode mixing 
resonance (|6aj) has been noted [lCj, 18 -23| but seems to 
have received less attention due to the fact that the main 
effect of interest has been particle creation. 

Mode mixing without particle creation is known in 
quantum optics as a passive transformation 27 1, imple- 
mented experimentally by passive optical elements such 
as beam splitters and phase plates. The oscillating cav- 
ity can hence theoretically be tuned to act as a beam 
splitter — a known quantum gate in continuous vari- 
able systems. Moreover, the entang ling power of passive 
transformations is well understood [l7, 28 "13 and is di- 
rectly applicable to our system. For example, the mode 
mixing generates entanglement from an initial Gaussian 
state only if this state is squeezed (25| |3C|. We empha- 
sise that while two-mode squeezing gates |l2 , 29|, 3(| and 
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other multipartite gates [3l| can be implemented by the 
particle creation resonance (|6b[) . mode- mixing gates are 
generated by the resonance l|6ap even when no particle 
creation is present. 

(3+1). — Let now be a real scalar field of mass 
fx > in a cavity in (3 + l)-dimensional Minkowski 
space, with Dirichlet conditions. The inertial cavity 
is a rectangular parallelepiped with fixed edge lengths 
L Xl L y and L z , and a standard basis of orthonormal 
field modes is indexed by triples (m,n,p) of positive in- 
tegers, such that the angular frequencies are u> mn p = 

\f V? + (nm/L x ) 2 + (irn/Ly) 2 + (irp/L z ) 2 . 

Acceleration in the cavity's three principal directions 

I 



can be treated as (1 + l)-dimcnsional, with the inert 
transverse quantum numbers just contributing to the 
effective mass. Acceleration of unrestricted magnitude 
and direction would require new input regarding how the 
shape of the cavity responds to such acceleration. To lin- 
ear order in the acceleration, however, boosts commute, 
and we can treat acceleration as a vector superposition 
of accelerations in the three principal directions in the 
cavity's instantaneous rest frame. 

To find the (3 + 1) Bogoliubov matrices, we first write 
the (1 + 1) matrices a and (3 in (|3b[) as a(M) and f3(M), 
indicating explicitly the dependence on the dimensionlcss 
mass parameter M ll|. We then define 



x^mnp,m' 'n' p' — &nn' $pp' Qmm f 



(i 2 L 2 x + w 2 n 2 {L x /L y f + n 2 p 2 (L x /L z ) 2 



zC^mnp.m'n'p' — $mm' $nn' Otpp' 



P > = 8 mm >8 pp , aw \ J H 2 L% + ir 2 m 2 (L y /L x ) 2 + n 2 p 2 (L y /L z ) 2 ^ , 



p?L\ + ^m 2 (L z /L x ) 2 + 7r 2 n 2 (L z /L v ) 2 



(7a) 
(7b) 
(7c) 



and we define analogously the matrices x /3, y f3 and z /3. 
Writing us = diag(w mnp ) and defining A and B as in (@|, 
equations §5§ generalise to 

A = iVa 3 (T)I J e -M-r-To) [ w a ] e fe(r-ro) ^ ( 8a) 
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-iui(T-To) 



, (8b) 



where [■,■] is the commutator, {•,•} is the anticommu- 
tator and (a x (r), a y (T), a z (r)) is the acceleration three- 
vector in the cavity's instantaneous rest frame. Once the 
acceleration is specified, A and B may be found from © 
by quadratures. 

For sinusoidal acceleration with angular frequency u r , 
the resonance condition is 



for A 



mnp,m n p 



for Hrnrtp,m' n' p' - — ^mnp ^m'n'p' j 



(9a) 
(9b) 



where in each case the difference in the quantum number 
in the direction of oscillation needs to be odd. 

Desktop experiment. — The mode mixing reso- 
nance angular frequency (|9ap can be significantly lower 
than the frequencies of the individual cavity modes. We 
outline an experimental scenario that optimises this low- 
ering. 

Setting // = 0, we assume that the quanta in the 
cavity have wavelength A <C min(L a; , L y , L z ) and have 
their momenta aligned close to the z-direction, so that 



(2/A) 2 w (p/L z ) 2 > (m/L x ) 2 + (n/Ly) 2 and uj mnp sa 
2tt/A + \-n\[(m/L x ) 2 + (n/L y ) 2 ]. We let the cavity un- 
dergo linear or circular harmonic oscillation orthogonal to 
the z-direction, with amplitude r x (r y ) in the x-direction 
(y-direction) . For motion in cc-direction, the mode mix- 
ing resonance angular frequency (|9a[) between modes m 
and mf , with m — m! odd, is 



\tt\L„. 2 \ m 2 



(m'Y 



(10) 



and it follows using (|5a|) . (|8a|) and the formulas for a(M) 
ll| that the mode mixing growth rate is 



^nmm / r x XL x ' 



(11) 



The lowest resonance occurs for m = 1 and ml = 2. Sim- 
ilar formulas ensue for the y-resonance, and for circular 
motion both resonances are present. 

As an experimental scenario, we first trap one or 
more quanta in the cavity, in modes whose momenta are 
aligned close to the z-direction. After a period of linear 
or circular oscillation perpendicular to the z-direction, 
a measurement on the quantum state of the cavity is 
performed. We assume that the resonance mode mixing 
dominates any effects due to the initial trapping and the 
final releasing of the quanta. 

We choose for example A = 600 nm and L x = L y = 
1 cm. The lowest resonance angular frequency is u> r f=s 
1.4 x 10 _2 m _1 w 4.2 x 10 6 s _1 , corresponding to an os- 
cillation frequency 0.7 MHz. With amplitude 1 ^m, (fTT]) 
and its y-counterpart gives jp|A rcs | «6x 10 2 s™ 1 , so that 
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the mode mixing coefficient grows to order unity within 
a millisecond. For longer periods the linear perturbation 
theory is no longer quantitatively reliable, but there is no 
reason to expect the mode mixing coefficients to decrease. 

For linear motion, oscillation of micron amplitude at 
megahertz frequency may be achievable by using ultra- 
sound to accelerate the cavity. Storing the quantum in 
the cavity for a millisecond could be challenging although 
recent achievements indicate that it may be feasible [32J . 

For circular motion, the threshold angular velocity 
iu r rj 4.2 x 10 6 s _1 w 4 x 10 7 rpm exceeds the angu- 
lar velocity 1.5 x 10 5 rpm achieved by medical ultracen- 
trifuges [33j . although by an amount that may possibly 
be bridged by a specifically designed system. 

Conclusions. — We have quantised a scalar field in 
a rectangular cavity that is accelerated arbitrarily in 
(3 + l)-dimensional Minkowski spacetime, in the limit 
of small accelerations but arbitrary velocities and travel 
times. The Bogoliubov coefficients were expressed as ex- 
plicit quadratures. For linear or circular periodic mo- 
tions, we identified a configuration in which the mode 
mixing resonance frequency is significantly below the fre- 
quencies of the cavity modes. If quantisation of the elec- 
tromagnetic field in an accelerated cavity is qualitatively 
similar, the mode mixing effects are within the reach of 
a desktop experiment that appears achievable with cur- 
rent technology in its mechanical aspects, if not yet in the 
storage capabilities required of a mechanically oscillating 
optical cavity. 

We also anticipate that the particle creation and mode 
mixing effects are not qualitatively sensitive to the de- 
tailed shape of the cavity, and this freedom could be 
utilised in the development of a concrete laboratory im- 
plementation. The experimental prospects could be fur- 
ther improved by filling the cavity with a medium that 
slows light down [34| . 

We underline that our experimental scenario does not 
involve significant particle creation. Nevertheless, it in- 
volves significant mode mixing. This mixing acts as a 
beam splitter quantum gate, creating or degrading entan- 
glement in situations where particles are initially present. 
We anticipate that observations of entanglement will gen- 
erally provide opportunities for experimental verification 
of the DCE that are complementary to observations of 
fluxes or particle numbers [§]. 
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